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BUCKLING TESTS VflTH A SPAR-RIB GRILL* 
By Josef tfeinhold 

SUMMARY 



The present report deals with a comparison of mathe- 
matically and experimentally defined "buckling loads of 
a spar-rib grill, on the assumption of constant spar sec- 
tion, and infinitely closely spaced ribs with rigidity 
symetrical to the grill center. The loads are applied 
as equal bending moments at both spar ends, as compres- 
sion in the line connecting the joints, and in the spar 
center line as the assumedly uniformly distributed spar 
weight* The evaluation formula is valid for a deflec- 
tion following a half wave, and confirmed by the tests 
within a scope sufficient for practical purposes. The 
formula is applicable with safety in the ratio bounded by 
by experiments between the end moments and the compres- 
sive force, or between the coefficients _s, , and m 2 » 
respectively. 

An effect of the finite aspect ratio of test spars 
(8 : 100) is practically not noticeable. Eor smaller 
compressive forces and greater end moments the coeffi- 
cients of the evaluation formula should be computed again 
for further deflection forms, in order to find further 
buckling conditions, which, of course, yield lower loads 
at the stability limit. For small spar weight a buckling 
condition exact for vanishing weight, which is written 
down for arbitrary wave number of the deflections, may be 
employed for the first comparison, which wave number 
gives the lowest buckling load. An extension of the en- 
ergy equation and consequently the achievement of a buck- 
ling condition for any degree of rib rigidity along the 
rib axis and other load distribution is readily accom- 
plished. 



*Uber Kippversuche mit einen Holn-Eippenros t . Luftfahrt- 
f or s Chung, vol. 17, no. 3, March SO, 1940, pp. 76-81. 
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Notation 






Mi. 


Ms* 


bonding r.onont In the spar en 


ds , nkg 




Si. 




compressive force in line connecting the 
kg 


joints , 


ft 

5 11 




spar weight, acting in bean axis, kg 




A 

A i> 


A 


rigidity of transverse bendin 






Ox, 


c 3 


torsional rigidity of spars, 


kgn 2 






Ar 


sun of rigidity of all ribs f 
plane of the grill, kgn s 


or buckling 


in the 




-D L 


sun of rigidity of all ribs for buckling 
angle to plane of grill, kgn 2 


at right 




A Ar 


strain energy of ribs for the spar length 
buckling in the plane of grill, kgn/m 


1 by 




-Br 


strain energy of ribs for the 
buckling at right angle to 

kgn/n 


spar length 1 by 
plane of grill, 


a U 


2 5 3 

il 


nondinens ional free value 
half distance of bean axes, n 








• 

.1 


distance of spar joints (spar 


length) , n 








end nonents of a rib for spar 


length 1 , k 


gn/n 




n 


half— wave nunber in the cours 


oof n • P i i 


and p 


° 1, 


c 2 


nunerical values of part integrals in the 
e quat ion 


energy 



x longitudinal coordinate, counted fron a joint, n 

y lateral deflection of bean axes at distance x 
fron joint, n 

*Subscripts l and 2 denote - if used alone - that the 
particular quantity belongs to spar 1 and-* spar 2, respective 
ly. Without subscript the quantities apply to both spars 
unless stated otherwise. 
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3-l, P 2 buck! i Jig angle (angle of rotation of spar sec- 

tions at right angle to "'plane of grill) 

cci , a 2 influence figures 

I = x/ 1 nonclimens ional longitudinal coordinate 



Ml + A 2 

T| = y/\ J «■» ncndinens ional lateral deflection 

v C 1 + C 2 

3:Io:idinensional quantities in the energy equation: 

2 Mi 1 v 2 M 2 1 

m i= : — . ra 



fUi + A 2 )(C 1 4- C 2 ) ' 7^ + A 2 )(C 1 + C 2 ) 

(s a + s 2 ) i 2 



6r 



/(Ai + A 2 )(Ci + C 2 ) ycXi + A 2 ) (0 : ! + C 2 ) 



2 Ar 1 2 Sr i 

D = 



(A 1 + A 2 )h 1 ~ ( C x + C 2 )h 

2 C, 2 C 2 

M»i = Jz — , M» 2 = 

C l +C 3 C l +C 2 



a x . - a 2 " - a 2 3 a 3 - a 4 



:::r?uODUCTiOff 



Earlier reports give buckling loads which were de- 
rived under specified simplifying premises and assumptions 
(references 1 and 2)-. In particular, the spars were as- 
sumed as snail rails with vanishing thickness and equal 
constant section for both spars. It is further assumed 
that the rib rigidity is constant and ' unif o rmly distrib- 
utable across the length of the spar. (Infinitely closely 
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spaced ri"bs as substitute, of f-inite number of ribs..) End 
"bending moments and end' compres sive forces were assumed 
as spar loads . 

The following describes an inves t i gat ion which enables 
a comparison of computed and experimentally defined' buck- 
ling loads of a spar-rib grill. The . exper Lmental grill 
has dissimilar spars which are not "very small" (8/100 mm), 
the number of ribs, 8, ■ i s c omparat i vely low). The rigid- 
ity of the ribs is no longer constant across the axis, 
although not symmetrically distribiited with resrject to the 
rib center. In the evaluation formula, which ordinarily 
is derived under the same general assumptions from the 
energy equation as before (reference 2), the dissimilarity 
of the spars, the cited course. of rib rigidity, and the 
spar weight are allowed for. 



Experimental Arrangement* 

The spars of the model grill are two pine boards 
8 by 100 by 3250 mm, selected for grain, absence of knots, 
etc., and the ribs (7 by 7 by. 1000 mm) of the same grade 
of pine with reinforcement strips glued at the ends (fig. 
2). Spar and rib end are fastened together by a 12-by-12- 
by-0.5-mn hard drawn brass angle and five 3-r.ir:: bolts* The 
spaci?ig of the spar axes from each other amounts to 2h ~ 
1025 mm. A hardened steel cylinder drilled out with a 
taper from each side and fitting tightly in the spar end 
reinforced by plywood shins, serves as nodal point of the 
spar ends. The annular edge within the cylinder has an 
3-mm diameter and rests against a hardened round steel bar 
of 7-min diameter. This in turn is fastened on one grill 
end to the support , on the other spar end the slidability, 
shown in figures 1 and 2, parallel with the line connect- 
ing the nodo-1 points, is provided. The end moment? are 
transmitted over fixed end levers built up of 10-by-l-mn 
flat brass and 10-by-10-by-l-mm angle brass, which project 
to 560-nm load application. For applying the axial com- 
pression, the fixed joint near the spars carries two ball- 
bearing pulleys over which the tension of the loading | 
weights is carried parallel with the connecting line of 
the joints to the movable joint. wooden ba,rs with metal 
sting fitted at the end levers of -spar 2 record the tan- 
gential inclination of the elastic line in the joints, on 
a little board with a sheet of paper. The length of the 
tracing levers from joint to metal sting is 1020 mm. The 
end moment existing on spar 1 as a result of the weight of 



*Seo fig. 1. 
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the end levers and of the overhanging , spar portion is 
0.082 mkg, and 0.100 mkg on spar 2 "because of the trac- 
ing levers. The effective weight on the joints (weight 
of spars plus weight of rits) is equalized for each spar 
at 1.21 kg* The "buckling angles in joints are preserved 
"by i-m wooden levers running parallel upward from the spar 
ends with a. loose pin or stud at the end parallel to the 
connecting line of the joints. The guide is adjustable 
perpendicularly to the lire connecting the joints. Fig- 
ure 4 shows the first set-up, where the spar plus the 
rib weight was upwardly compensated "by pulleys. But this 
compensation had to "be abandoned in favor of a more com- 
plicated evaluation formula "because the oscillations of 
the grill could not ho removed quickly enough after load 
changes. Soither figure 4 nor figure 3 represents the 
final method. * 



Evaluation Formula 

The energy equation for elastically identical spars 
with a load through U% 9 M 2 and Sj and S 2 and ribs 
of constant section (reference 2) reads: 



j = p |(M lPl + M 2 P ,) 7 - -(^-^ - 3 £) y .- 



I 

"» v 2 Ih 

* — (Pi 2 + PiP 3 •:- £ 2 2 ) + | C(p x ' 2 + 3 2 «2) + A y»s \&x 
Ih 2 J ( 

or in nondimens ional coordinates and coefficients 
1 

Jo 

+ i Oi ,a + * n ,,s l *t (2) 

2 J ... 

M x and M 2 are the principal moments "bending at point x 



*The tests were made in the German Technical Institute at 
Brunn, hut had to be interrupted during the general mobil- 
ization in September 19 38. The means for the tests had 
"been furnished by the National Research Council of Czecho- 
slovakia. 
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Accordingly, bjf virtue of the annexed loading schedule, 
VL 1 and Mg must "bo replaced "by, respect ively , 

Mi - —i-(3ct - x 2 ) and M 2 - -£(xl - x 2 ) 
21 2t 

For the derivation of the proportions with Ar and 
Br the "buckling of a riq at right angle to the plane of 
the grill due to the reaction j>arts transmitted at the 
joints, is expressed coiif crmally to figure 7, with 

P = <x x - m 3 a s ) , p 3 « M.(m 2 ai - »j co 2 ) (3) 

3r Br 

(The rib rigidity is ■— , the rigidity for spar length 1.) 
Then the energy content of the rih is 

A Br = Ji_ (^(p^ + p 2 2) + 2 a 3 P 1 P 2 ) (4) 

2 h Ua^ - a 2 '*') 

For buckling in the plane of the grill, it affords in the 

same way 

Ar , v 

AAr = TTT? r y' 2 (D) 

h I k<x 3 - a 4 ; 

if y 1 (reference 2) replaces p x and g 2 , and a 3 and 
a 4 replace ex x and a 2 . 

The sum then is 
I 



J = 



'0 



- - »*)) »* + («• - |f (»» - *•>) p 3 } y" 

2 hl(a 3 - a 4 ) J 2 h l( ai s - a 2 2 ) 

+ | C^* 8 + | 0 3 .6 2 ' 2 + J Uj + A 2 )y« 3 d :c (6) 
2 & & 

Ol 3 

Extension with , followed hy introduction of the 

Ox * G 2 

nondinens ional coordinates and abbreviations , gives: 
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The limiting conditions adapted to the described test ar- 
rangement read (case 1, reference 2): 

n(o) - n(i) = n"(o) - n"U) = o \ 

( (8) 
0 X (O) - = 3 3 (0) ■ Ml) - 0 J 

These are experimentally much easier to keep than others 
and give at the same time lower 'buckling loads in ratio 
to built-in spars. 

The related differential equations following from the 
disappearance of the first variation of equation (?) are t 

1 d s f 

2 d i 



pg - siU - l*))$ x + tm - frit - i 2 )) p 3 } 



+ n IV + (s. - La))! 11 = 0 (9) 



(as - GrMi - I 2 ))n" - n 3 # 2 « + i(Sa^+ Ms) = 0 f 



(10) 



After two integrations of equation (9) the integra- 
tion constants disappear as a result of the chosen limit* 
ing conditions, leaving the simpler relation 

i (mi - £iCt - i z ))px + J(a a - e 2 (£ -e»)) 0a 

-i- n M + (s. - £ 3 &)f| = 0 (9a) 



For vanishing spar weight equations (9) and (10) are satis* 
fied "by 



B X = 



Ki s i n n TT |i (5 2 s 



K 2 s i n n it | , 



n = K 3 ;S in n tt I (11) 



8 



NACA Technical Memorandum No* 950 



The final determinant remaining after the insertion of 
the solutions (ll) gives the "buckling condition 

2 ( rnr ) 2 1 ^ ( mr ) b i 

a, as f v /A ^ t«& , \ 
- 7^T5 t + v~ — dr ~ 1 J 

(mr)* + 2 £ x £ + (t x 2 - iV = 0 (12) 



The exact solution of the complete differential 
equations was foregone for reasons of the paper work in- 
volved, and the evaluation formula therefore derived hy 
presentation of suitably chosen approximating theorems in 
the energy equation. Chosen as such were: 



i. pi = ~a - e 3 ). p 2 = ~(g - e 3 ), n« - — ct - i s ) (13) 

1 2 3 

ii. 3i «iHe-2£ 3 +e 4 ). Pa=^!^| 3 +| 4 ), n"«4-(!-3| 3 +| 4 ) (14) 



In order to clarify the mutual relations of the coeffi- 
cients in the "buckling equation the numerical values of 
the part integrals in the energy equation are generalized 



'o 




1 




(17) 



J 

0 
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• : i " a , 

TfMt- ' 2)d 8 = "777" / *l B Mt- e 8 )dt=c 4 (18) 



B X R 3 J K 2 K 3 

o . si 



With these data the energy equation (7) supplies as buck* 
1 ins condition 



SUM til) [f - 2-^m l£ 



-I- £*, 



9 C4 



- 1 £3 — BLi&a ~ c7 fax«i. + BsjSs) + gigs 



C A 



C i c 



2 



* Ifu- i^)- 1 ) i^pvE- * '2 ( i x 8 - u 3 )i 2 £j } * o 



f (19) 



The numerical values of c lf C a • . • and their com- 
binations are given for approximations I and II in tho 
second and third lines of table I. 

A comparison of equation (19) with the exact solu- 



tion (.12)' discloses, then, 



and 



C3 



for tr 2 . Moreover, 



the limiting conditions considered here give for 
limoshonko f s (reference §) exact buckling factor 



Gr t 



>8. 3 



(20) 



for 



c , c 



the actual value i 28,3 



C4 



200. 22 3 . The third 



approximation and final evaluation formula employed is 

that in which §1* tl t and ii£ft are replaced by the" 
u 1 c 3 c 4 c 

exact values and |£ by the value 44.4535, the latter 



following fro m 



TT ' 



and 



c i c ■ 



C 4 < 



== 200. 22 3 and which 



probably differs very little and none s c ent ially from the 
true value 
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i ♦ t,l) ! !f- - ftft ♦ 2C 



m 
lTt : 



+ I (u^TT 2 + ^1) 



IT 



22.23 200.22 



&i&a 



TT 



TT 2 



44.45 



ft, £ 



200.22 



- l) { u-^tt 2 + 2£ 2 b + - U S ) ^J) = 0 (21) 



A-x • Oi, A 2 , C 3 were defined from buckling tests with one 
single spar each without ribs attached. The related evalu- 
ation formula follows from the above after insertion of 
Mj, = M 3 = M, Si a S 2 = S, a = 0, b = 0, and the expres- 
sions for _s, m 1 , and m 2 in the form 

iLJLi - B L? M i + ^ i 4 + c s i 3 = A c (22) 

tt2. 44.45 , 800.39 rr^ . 



Tost Procedure 

The first tests were made without the weights sus- 
pended from the end lever, one horizontal pull cord being 
used to transmit the principal part of the load, and then 
the other. After various attempts the locations of the 
cylinders slipped in the spar ends were finally established 
for which the deflections by eq\ial load near the stability 
limit were lowest. The inclination of the elastic line 
at the spar ends served, as already indicated, as deflec- 
tion component, the board with the sheet of paper being- 
pressed against the tracing lever after the set-up ceases 
to swing and the load is written alongside the tracing 
mark. The adjusting levers for end-buckling angle remained 
loosely in the guides. ?or further tests with this set- 
ting the grill was fitted with 7 kg end-lever loads on 
spar 1, on spars 1 and 2, and then on spar 2. The applied 
load at the end levers could be shifted transversely to 
the beam axis, and so far that after application of the end- 
lever load the resetting lever over the joints was* loose 
again in the guides. Thus the end load was applied in 
consistently smaller stages till the lateral deflection in 
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spar center reached about 100 mm. Greater deflections 
produced appreciable permanent deformations. Up to 100 
mm lateral deflection in spar center the permanent de- 
flection, after applied end load, was at the most 1*5 mm 
at the pointer tip. For the same reasons the recording 
of a supercritical "branch of the stress-strain diagram, 
as by Prandtl (reference 4) was omitted. The height of 
the critical end load had to be determined, therefore, 
as hyperbolic asymptote height by known asymptote direc- 
tions and three given hyperbola points (reference 5). 
The force deflection diagrams taken at both ends of spar 
2 differ very little, which is indicative of the symmet- 
rical course of the deflections. Tests with 10 -kg end- 
lever load left perceptible permanent deformations, es- 
pecially in the spar-rib connections and the symmetrical 
course of the deflections disappeared. The evaluation 
formula would then give too high comparative values. Ex- 
periments and comparisons in this direction with an im- 
proved set-up are under way. 

Table II contains the end -moments Mi , M 2 applied 
at the grill, and the related end compressive forces de- 
termined from the deflection plots. Figure 8 is the dia- 
gram showing the comparatively greater departure from the 
presumed symmetry of deflection, whose equilibrium curves 
have the relatively greater distance from the asymptotes. 
Sj + S 2 =5 11.26 kg is the sum of the "scale loading" of 
10.8 kg determined by the asymptote construction and the 
"scale weight" of 0.457 kg identical for all the tests. 

Spar Rigidity A ly G 1% A 2 , C 2 

Following the buckling tests with the model grill 
the ribs were removed and each spar tested at three dif- 
ferent end moments. The effective weight of each spar was 
G.980 kg. The three load systems are appended in table 
III. Equation (22) gives for A x , Oi, A 2 , C 2 three equa- 
tions each which, interconnected, give each three pairs 
of rigidity A and C. Their averages 

A 1 = 4.961 m g kg Oi = 0.605 m s kg 

A a = 4.978 efjgg C 2 = 0. 655 m^k g 



served as a basis in the subsequent calculations. 
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Determination of Rigidity Factors ^2^L> and - ^3^1 



Two each of the eight ribs spaced 15 cm apart were 
attached to a 20-cm long piece of spar in the same way as 
on the grill and supported and loaded. The torsion anglei 
over the joints were determined by means of mirror read- 
ing and compiled in table IV. It shows the extent to 
which the symmetry assumption is satisfied for the course 
of rib rigidity along the rib axis. Equation (3) and the 
equations defining b, and £ s give 



Hi ■ 



2 mj_ \ s ^ 1 



(C x + C 3 ) (Px 3 - 3 2 s ) 



and h£ 3 = 



2 m 1 l 3 3 2 



(Ci + C 2 ) 



P 2 



8 ) 



(23) 



and in our case for the mean values of Pi, p 2f and £ 2 * 

t x b = 50. 6 , £ 2 b = 2 6.5 



I = 0.0423, 1 = 0.0218 



Test 4 without end-lever load affords ts&i which, 
gives Sj, + S s = 16.52 kg. The evaluation formula (21) 
then ^ives 



Q a a 5.14 



According to the definition of £ 3 a, ^jjb and £ s b_ should 
be - equal rigidity of ribs and their attachments in and 
at right angles to the plane of the grid being assumed - 



S3- b»a Ai + As S.949 

The much lover actual value of 5.14 is due to the fact 
that the rib fittings are much more elastic with respect 
to buckling in the plane of the grill than at right angles 
to the plane of the grill. It might be noted in passing 
that the eight ribs had boon selected from about 80 others. 
The sequence of arrangement on the spar is that of table 
17, hence the pair closest to the average value is situat- 
ed at spar center where the greatest twist occurred. 
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Comparison of Theoretical and Experimental Values 

Table II gives the end moments, the critical compres- 
sive forces found from the force-deflection diagrams, the 
mathematical compressive forces, and the difference "be- 
tween the compressive forces in percent. In view of the 
uncertainty affecting the determination of the critical 
compression forces from force-deflection diagrams .and par- 
ticularly "because of the dissimilarity of the ribs, the 
accord "between theory and test may "be termed sat i sf ac t ory • 
Since the evaluation formula for the snail rail is exactly 
valid, the critical experimental compressive forces should, 
on the whole | exceed the theoretical, which is the case 
in test 1 only, For the present 8:100 aspect ratio a pure 
moment loading through M x and M 2 • with and without al- 
lowance for finite aspect ratio would give a difference 
of ah out 1 percent, which, however, is covered "by the un- 
certainty of the coordination. In test 3 the difference 
between the theoretical and the experimental critical com- 
pressive force with -6»2 percent is greatest compared with 
the two previous tests. This can "be associated with the 
previously stressed maximum departure from the symmetry 
hypothesis for the deformation and this in turn with the 
dissimilarity of the ribs. 



Translation "by J. Vanier, 
Hfat ional Advi s ory Commit tee 
for Aeronautics . 
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TABLE I 



Approx- 
una oion 


c i 


^2 


I ' " 
c 

°3 


c 4 


— , 

it 


1 


c ! c 2 
c 4 


Cl 

c 3 


I 


1 


1 


1 

17 


1 


10 


46.6 

I 


217.7 


9.86235. . 




30 


3 


35x144 


140 


II 


31 


17 


691 


37 


9.870P68. . 


45.486. . 


209.607. . 


9.86975. . 


13x85 


35 


09x1400 


3465 


III 










TT 2 


| 

44.4535 


200.223 
j . . 


TT 2 



TABLE II 



Critical Outside Loads of Grill 



Test 


Mi 


K« 


(-i + s 2 ) test j 


1 ° ' theory 


P * 




(mkg) 


(mkg) 


(kg) 


(kg) 


(percent) 


1 


0.082 


4.02 


... 

12. 29 


11.95 


2.8 


2 


4.002 


4.02 


11.56 


11.82 


-2.2 


3 


4.002 


.10 


11.26 


12.01 


-6.2 


4 


.032 


.10 


16. 52 

1 


16.52 





■^Difference "between experimental and theoretical value in 
percent of computed value. 



r 
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TABLE III 
Buckling Loads for the l n roc Spars 



Spar 1 


Spar 2 


M 

(mkg) 


^ S 
(kg) 


M 

(mkg) 


S 

(leg) 


0.10 


5. 35 


0.10 


5.38 


1.21 


4.14 


1 . 21 


4.22 


1.78 


1.95 


1.78 


2.115 



TABLE IV 

Rib Torsion Angle for mi 4" and m s 4* = 0.0813 ink 

4 4 



Sit 


Pi 


Pz 


B 1 


fa 


1 

2 


0. 0435 


-0 .0225 


0. 0228 


-0.0438 


3 

4 


. 0431 


- .0222 


.0218 


- .0427 


5 


. 0419 








6 


- .0214 


.0215 


- .0416 


7 










8 


. 0409 


- .0212 


.0210 


- .0409 


Average 


. 04235 


- .021825 


.021775 


- .04321 
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Pigs. 1,2, 3. 4, 8. 




Mil 



± jrM Spar 2 



Figure 1.- Experimental set up. 





Figure 4 



Figure 2.- Spar support, rib fitting and 

adjustment of end buckling angle, 



First experimental 
6et up. 




Figure 3.- Pulley at fixed joint. 



s^Sf^kg Scale weight. 

no 



10,0 





Mb 


















Test No 

ff t - 002 mki, 
H z = 0,10 - 
S f tS 2 "11,25 


.3 

r 

7*} 






o y'(0). 1020 m 
•y'(l). * ■ 





Figure 8.- Load-deflection 
diagram for test 
No. 3. 
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Jigs. 5 




figure 6.- Bpar-bending loads. 
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Figs. 7,9. 




Figure 9.- Test rigging for determining B r , 



